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Gauge invariant investigation of the nature of Confinement 
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We observe a strong correlation between the decrease in the number of action density peaks in SU(2) Yang- 
Mills configurations with cooling and that of the string tension. The nature and distribution of these peaks is 
investigated. The relationship with monopole currents after the abelian projection is also considered. 



1. Introduction 

In this paper we report on some of our more 
interesting results within our program to analyse 
the physical content of the SU(2) Yang-Mills vac- 
uum. Our final goal is to achieve a comprehen- 
sive description of the structure of the vacuum 
and, hence, an understanding of its more salient 
properties: Confinement, chiral symmetry break- 
ing, etc. A similar goal is pursued by other groups 
with contributions collected in these Proceedings. 
Our work, nonetheless, focuses on gauge invariant 
structures and is hence different, but not, neces- 
sarily antagonistic, from attempts based on the 
abelian projection. 

Our methodology is, in broad terms, standard. 
We have Monte Carlo generated 56 configura- 
tions on an 8 3 x 64 lattice with twisted boundary 
conditions in space ( m = (1, 1, 1)) and Wilson 
action with [3 = 2.325. Then, we have cooled 
these configurations using over-improved cooling 
[1] and e = —0.3. Finally, we have studied the 
gauge-invariant properties of the resulting config- 
urations. Before proceeding to describe our re- 
sults, let us clarify some technical points about 
our method. 

The spatial size of our lattice in physical units 
is l s = 1.25 fm. ( with a = 5 fm~ 2 ). Although at 
this size there are probably non-negligible finite- 
size effects, we claim that qualitatively the struc- 
ture of the vacuum is essentially the same as for 
infinite volume. We base this belief on the results 
of our previous M.C. simulations at this same 
point [2], where we measured the string tension a 
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from correlations of Polyakov loops with different 
windings in space. We got results consistent with 
the infinite volume string tension with up to 10% 
differences and errors. Both the l s dependence 
and the equality of the estimates from correla- 
tors with different windings (electric fluxes) are 
in agreement with the string formation picture. 

Another point of concern is the effect of cooling. 
Contrary to other authors, we do not try to argue 
that the results obtained are unaffected by this 
technique. Rather, we have carried an indepen- 
dent study to understand the effect that cooling 
might have on the results. We have looked at pre- 
viously prepared configurations, cooled them for 
a while, and observed the outcome. Since cooling 
is a local minimization technique it is clear that 
classical solutions will sooner or later dominate 
the configurations. Thus, we first examined the 
effect on individual instantons. Contrary to the 
case of the continuum, instantons on the lattice 
are not invariant under cooling. The degeneracy 
of zero-modes is broken by artifacts and cooling 
produces a motion along this valley. In the stan- 
dard cooling algorithm, instantons become small 
and eventually disappear through the holes of the 
lattice. Our proposal of over-improved cooling 
makes use of the fact that one can change the size 
and the sign of artifacts to make this effect smaller 
and eliminate the instability of sufficiently large 
instantons. In Ref. [3] we advocated that this 
variant cooling method could be used to obtain a 
better determination of the topological suscepti- 
bility. (One poster in this Conference [4] presents 
a more refined proposal within this idea.) In this 
work we have chosen e = —0.3. With this choice, 
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Figure 1. The topological susceptibility as a func- 
tion of cooling step for e = -0.3 (closed squares) Figure 2. Distribution of the peaks in the p x (ver- 
and ordinary cooling (white circles). tical) and p 2 (horizontal) plane. 



large instantons are stable, but at the same time 
motion along the zero-mode valley is not too fast. 
Concerning individual instantons, hence, we con- 
clude that our cooling essentially simply removes 
those which are smaller than 2 lattice spacings. 
In Fig. 1 we show the topological susceptibility 
as a function of the number of cooling steps for 
e = —0.3 and for naive cooling (e = 1). 

We have also investigated the effect of cooling 
on instanton-antiinstanton pairs. We have ob- 
served that cooling produces the annihilation of 
such a pair, although the rate at which this hap- 
pens depends widely on separation and sizes of 
the instantons. This mechanism is independent 
on the value of e and, hence, cannot be elimi- 
nated. We should then understand its effect and 
correct for it. 

Now we should consider which are going to be 
our main observables. We will look for peaks (i.e. 
local maxima) in the action density, its electric 
and magnetic parts separately and the topologi- 
cal charge density. For each peak we will record as 
attributes, its location, its height and its width, 
obtained from the first neighbours of the peak po- 
sition. We aim at an understanding of the nature 
of these peaks and a description of their distri- 
bution from our data. Some of these peaks will 



be instantons, of course. Since they depend on a 
single scale p, for them the height and the width 
are related. On the continuum we can form two 
independent functions of height and width p\ and 
P2 as follows: 
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where S(x) is the action density for the configu- 
ration with a peak located at x = and AS'(O) 
its 4-dimensional laplacian at this point. For a 
BPST instanton p\ and p^ give the same value 
and equal to the instanton size parameter p. 

2. Results 

In a previous publication [3], we showed the be- 
haviour of the string tension and the topological 
susceptibility as a function of the number of cool- 
ing steps. As reported by many authors before, 
both quantities decrease under the naive cooling 
method. As mentioned previously, if we use over- 
improved cooling and e = —0.3 the decrease in the 
topological susceptibility is eliminated (see Fig. 
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Figure 3. The quantity A'(5) = (open circles) 
versus the number of cooling steps. The effective 
string tension in fm~ 2 units (closed circles) 
is also shown. 



Figure 4. We plot the average action of a config- 
uration per peak divided by 4ir 2 for each cooling 
step. 



1), but however that is not the case for the string 
tension, which at 50 cooling steps has dropped to 
40% of its uncooled value. Let us clarify what is 
meant by string tension for the cooled configura- 
tions. We may measure an effective string ten- 
sion a(t) through correlations of Polyakov loops 
at distances t and t — I. This quantity shows first 
a plateau for large but finite t, which is the quan- 
tity we are referring to. Nonetheless, if t is large 
enough (t > N c ), as argued by Teper, we should 
recover the uncooled string tension. 

In the present work, as mentioned previously, 
we studied the nature and distribution of peaks. 
The first result concerns the degree of self-duality 
observed. The peaks in the electric action den- 
sity coincide in 90% of the cases with those in 
the magnetic density (or the topological charge 
density). The coincidence is in number, location 
and height. Usually, if spurious or non-symmetric 
peaks appear, they seem to be associated with 
transient processes of instanton-antiinstanton an- 
nihilation. Hence, from now on we will not specify 
which peaks we are referring to, since they all give 
the same answers. Our next observation, is that 
the number of peaks decreases with the number 
of cooling steps N c . From 25 to 50 coolings this 



number has dropped by a factor of 2. It is quite 
natural to interpret this decrease as due to our 
previously mentioned instanton-antiinstanton an- 
nihilation mechanism. But, are our peaks BPST 
instantons? To see this, we checked the equality 
of p\ and p2 for our peaks. In Fig 2. the con- 
tour plot in the P2P1 plane is given. Most peaks 
are lying along the diagonal, with deviations due 
to finite-lattice corrections. This is a non trivial 
fact, which need not be the case even for self-dual 
peaks. We conclude that our peaks are locally 
instanton-like. They need not be BPST instan- 
tons, since there are other configurations showing 
the same behaviour. In particular, the Q = ^ in- 
stanton of Ref. [5] is known to be present and to 
behave in this way. 

An important observation that we made, is the 
existence of a complete correlation between the 
decrease in the density D of peaks and that of 
the string tension. We may form a dimension- 
less quantity K = and study its evolution 
with the number of cooling steps. The result is 
invariant under cooling to a 2% level from 20 to 
50 cooling steps. Being more precise, we may 
write K(t) for the quantity obtained from the ef- 
fective mass &(t). Our statement applies for all 
values oft from 2 to 6. In Fig. 3 we show A'(5). 
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The value of in fm~ 2 units is also plotted, 
to show how the string tension decreases in this 
same interval. If we fit a constant value num- 
ber to the values obtained from 20 to 50 cooling 
steps, we get A(5) = 2.18 ± 0.205 ± 0.027 where 
the first error is statistical (for a single cooling 
step), and the second is the maximum difference 
observed in absolute value between the constant 
and the value for all cooling steps from 20 to 50. 
For other t we get: A(2) = 0.939 ± 0.024 ± 0.027 
, A(3) = 1.463 ± 0.048 ± 0.037 , A(4) = 1.881 ± 
0.097 ± 0.036, A(6) = 2.331 ± 0.427 ± 0.028, 
A(7) = 2.236 ± 0.813 ± 0.048. Notice that for 
the larger values of t, K(t) tends to a constant 
(which is what we call string tension). 

In Ref. [6,2] we proposed an scenario for the 
Yang-Mills vacuum based on a gas of self-dual 
\Q\ = i lumps with mean separation d around 
0.7 fm. This picture predicts precisely that the 
string tension should scale with the square root 
of the density d~ 2 . Using a = 5/m -2 and 
d = 0.7 fm, we get K = 2.5, which is close to 
the estimates from our data. Our model has obvi- 
ously a crucial test: the observed lumps should be 
associated with ±-| topological charge. To check 
this point, we computed the total action of each 
of our configurations and divided it by the num- 
ber of peaks of that configuration. The value ob- 
tained was in all cases close to Ait 2 (the value for a 
self-dual |Q| = ^ object). Averaging over all con- 
figurations, we obtain the mean action per peak 
as a function of cooling steps, which is plotted 
in Fig. 4. Notice the small errors, reflecting the 
small dispersion of the data. The value is also 
fairly stable under cooling given the large varia- 
tion of numerator and denominator in the 20 to 
50 step region. We lack a simple algorithm that 
would tell us what peak is associated to Q = ^ 
and which to a normal instanton. Nonetheless 
we possess a method which is non-local and quite 
lengthy to be used, but can be used for a few con- 
figurations. These control configurations at 100- 
200 cooling steps are very clearly seen to be made 
of a majority of |Q| = ^ lumps and a few integer 
charge ones. For fewer number of steps the iden- 
tification is less beautiful and requires taking into 
account transient instanton-antiinstanton annihi- 
lation states as well. 



We conclude with a summary and comments. 
Our data provides clear evidence for the corre- 
lation between the string tension, as measured 
through Polyakov loops, and the density of peaks. 
The latter decreases with cooling via pair annihi- 
lation (also |Q| = i pairs annihilate), and drives 
the former. The observed correlation is to be ex- 
pected within our picture of the vacuum based on 
a gas of |Q| = ^ lumps. Some evidence for the 
presence of these lumps has been given, both in 
an statistical fashion as by monitoring individual 
configurations. A more detailed description of our 
results can be found in Ref [7]. Let us mention 
to conclude, that we have looked (together with 
M. Polikarpov) at what results from the |Q| = -| 
configurations after the abelian projection. When 
using Maximal abelian gauge, we found a single 
monopole line going through the center of the 
lump. This is again the case if one takes ex- 
tended monopoles of size 2a. For other gauge 
fixings (Polyakov lines or plaquette) a more in- 
volved structure of monopole lines is found. More 
work in this direction is in progress. We conclude 
that, presumably, a gauge-invariant description 
of Confinement is not in contradiction with the 
Dual Superconductor picture. 
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